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ABSTRACT 

There are now many known exoplanets with Msini within a factor of two of Neptune's, 
including the transiting planets GJ436b and HAT-P-llb. Planets in this mass-range are 
different from their more massive cousins in several ways that are relevant to their radiative 
properties and thermal structures. By analogy with Neptune and Uranus, they are likely 
to have metal abundances that are an order of magnitude or more greater than those of 
larger, more massive planets. This increases their opacity, decreases Rayleigh scattering, 
and changes their equation of state. Furthermore, their smaller radii mean that fluxes from 
these planets are roughly an order of magnitude lower than those of otherwise identical 
gas giant planets. Here, we compute a range of plausible radiative equilibrium models of 
GJ436b and HAT-P-llb. In addition, we explore the dependence of generic Neptune-mass 
planets on a range of physical properties, including their distance from their host stars, their 
metallicity, the spectral type of their stars, the redistribution of heat in their atmospheres, 
and the possible presence of additional optical opacity in their upper atmospheres. 

Subject headings: equation of state - line: profiles - planetary systems - radiative transfer 
- stars: individual GJ436, HAT-P-11 - astrochemistry 



1. INTRODUCTION 

Although in the early days of exoplanet research 
large, massive planets (gas giants) were discovered 
in disproportionate numbers, lower mass objects, in- 
cluding super Earths (1 to '^10 times the mass of the 
Earth) and ice-giants, are now being found with in- 
creasing frequency.^ The latter population, consist- 
ing of "Neptune-mass" bodies (those with masses 
within a factor of a few of Neptune's), currently 
numbers several dozen, and is growing rapidly. 

Exoplanet transits, first observed by Henry et al. 
(2000) and Charbonneau et al. (2000), allow precise 
measurements both of planet masses (by breaking 
the degeneracy between mass and inclination angle) 
and of radii. These two crucial pieces of data pro- 
vide constraints on models of structure, evolution, 
and bulk composition (Ibgui & Burrows 2009; Bur- 
rows et al. 2007a, 2000; Guillot 2008, 2005; Guillot 
& Showman 2002; Baraffe et al. 2005). As a result, 
the best-studied among the Neptune-mass popula- 
tion are the transiting planets GJ436b (1.33 times 
the mass of Neptune, and 1.27 times its radius; Cac- 
eres et al. 2009; Bean & Seifahrt 2008; Bean et al. 
2008; Torres 2007; Deming et al. 2007a; Maness 

Electronic address: dsp@astro.princeton.edu, bur- 
rows@astro.princeton.edu, ibgui@astro.princeton.edu, 
hubeny@as.arizona.edu, milsom@physics. arizona.edu 
^ See the catalog at http://exoplanet.eu. 



et al. 2007) and HAT-P-llb (1.507\fNcp, 1.30i?Nop; 
Bakos et al. 2009; Dittman et al. 2009). 

Observations of an exoplanet can constrain the- 
oretical models. The marriage of data with theory, 
therefore, allows us to learn about the physical con- 
ditions on these distant worlds. For example, obser- 
vations during transit can constrain various prop- 
erties, including both the atmospheric composition 
and structure at the terminator (Seager & Sasselov 
2000; Hubbard et al. 2001; Charbonneau et al. 2002; 
Desert et al. 2008; Redfield et al. 2008) and the mo- 
tion of the atmosphere (Brown 2001; Spiegel et al. 
2007). This paper, however, focuses on the diagnos- 
tics of planetary structure that are available from 
observations of a planet's disk, not of its limb. The 
disks of those planets whose orbits are aligned pre- 
cisely enough that they transit may be studied in 
exquisite detail. Half an orbit after transit, they 
undergo secondary eclipse, when they pass behind 
their host stars. Immediately before and after sec- 
ondary eclipse, a planet is in its "full-moon" phase, 
when its reflected and emitted light are at their 
maximum from our vantage. There have been a 
variety of ground-based and space-based secondary- 
eclipse observations of exoplanets in the last several 
years (Richardson et al. 2003; Snellen 2005; Dem- 
ing et al. 2007a, b; Charbonneau et al. 2008; Knut- 
son et al. 2008, 2009; Alonso et al. 2009b,a; Sing 
& Lopez-Morales 2009). In the last decade, several 
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groups have produced a progression of theoretical 
models of the atmospheres of Jupiter-mass exoplan- 
ets (Sudarsky et al. 2000, 2003; Hubeny et al. 2003; 
Burrows et al. 2005; Fortney et al. 2006; Burrows 
et al. 2008a; Fortney et al. 2008; Burrows et al. 
2008b; Spiegel et al. 2009; Showman et al. 2009), 
including calculations of both optical albedos and 
infrared emergent spectra. These calculations are 
now being used in interpretative studies of the sec- 
ondary eclipse observations.^ 

A critical difference between the gas giant and 
the ice-giant planets motivates the present investi- 
gation of theoretical models of lower mass objects. 
The gas giants of our solar system - Jupiter and 
Saturn - have metal abundances not more than 
a few times solar (Matter et al. 2009; Saumon & 
Guillot 2004). However, Neptune, Uranus, and, 
presumably, extrasolar planets in the same mass- 
range, are thought to have bulk metallicities of 30- 
45 times solar (Figueira et al. 2009; Matter et al. 
2009; Guillot & Gautier 2007). Increased mctal- 
licity leads to increased opacity. Reduced relative 
hydrogen abundance might also reduce the impor- 
tance of Rayleigh scattering in the atmospheres of 
Neptune-mass planets. 

In the last several years, infrared observations 

by the Spitzer Space Telescope, in conjunction with 
theoretical models by our group and others, have 
suggested that several giant planets have thermal 
inversions, wherein, above a relative minimum, the 
atmospheric temperature increases with altitude. 
Hubeny et al. (2003) suggested that, if there is an 
additional source of upper atmosphere optical opac- 
ity, incident stellar irradiation could lead to pre- 
cisely this type of atmospheric structure. The plan- 
ets whose emergent spectra have suggested ther- 
mal inversions include HD 209458b (Burrows et al. 
2007b; Knutson et al. 2008), HD 149026b (Fortney 
et al. 2006; Burrows et al. 2008a), TrES-4 (Knut- 
son et al. 2009), XO-lb (Machalek et al. 2008), 
X0-2b (Machalek et al. 2009), X0-3b (McCullough 
et al. 2008), and, perhaps, v Andromeda b (Bur- 
rows et al. 2008a). Various authors have suggested 
that titanium oxide (TiO), as a strong optical ab- 
sorber, might provide the extra opacity that is 
needed to produce the inferred inversions (Fortney 
et al. 2008; Showman et al. 2009). However, Spiegel 
et al. (2009) argue that, without extremely vigor- 
ous macroscopic mixing, a heavy molecule such as 
TiO will settle to a level that is too deep for it to 
contribute to an inversion. Furthermore, Spiegel 
et al. (2009) show that several planets for which 
inversions have been inferred ought to have cold- 
trap regions deeper in their atmospheres that are 
cool enough for titanium to condense and rain out, 
therefore requiring even more vigorous mixing if 
TiO is to survive up to the millibar levels where 

^ Comparable calculations have been made in the context 
of Super-Eaxth planets (Miller-Ricci et al. 2009; Kaltenegger 
et al. 2009). 



it would be needed.'^ Zahnle et al. (2009) suggest 
that sulfur photochemistry provides another avenue 
for achieving the additional upper-atmosphere op- 
tical opacity that is needed to produce inversions. 
The ultimate cause of inversions remains unknown, 
and might differ from one planet to another. Nev- 
ertheless, sources of extra optical opacity might be 
related to metal abundance, and so it is reasonable 
to wonder whether some highly irradiated Neptune- 
mass planets might have thermal inversions. 

In order to survey the range of plausible struc- 
tures of Neptune-mass exoplanets, we compute a 
variety of model atmospheres. We explore how 

the emergent infrared fluxes and optical albedos of 
GJ436b and HAT-P-llb depend both on the redis- 
tribution of heat in their atmospheres and on the 
possible presence of an extra source of atmospheric 
opacity (following the treatment in Burrows et al. 
2008a and Spiegel et al. 2009). We also examine 
how "ice-giant" atmospheres, and their optical and 
infrared fluxes, depend on their distances from their 
host stars, their atmospheric metal abundances, and 
the spectral types of their hosts. By "ice-giant," 
we are simply referring to a planet with approx- 
imately Neptune's mass and radius; whether the 
atmosphere is bounded below by ices is irrelevant 
to our analysis. Otir models are of planets with 
hydrogen/helium-dominatcd atmospheres that are 
metal-enriched relative to solar by a large factor. 
We note that some planets currently classified as 
"Super-Earths" might share this structure. 

The remainder of this paper is structured as fol- 
lows: In §2, we present our numerical techniques for 
computing the atmospheres of "ice-giant" planets. 
In §3, we describe the different models that we ran. 
Section 4 contains the results of our calculations. 
Finally, in §5, we summarize our findings. 

2. NUMERICAL METHODS 

We calculate radiative equilibrium models of ir- 
radiated planetary atmospheres models. As in our 
other recent studies, we use the code COOLTLUSTY 

(Hubeny et al. 2003; Sudarsky et al. 2003; Burrows 
et al. 2006, 2008a; Spiegel et al. 2009). This code 
is an variant of the code TLUSTY (Hubeny 1988; 
Hubeny & Lanz 1995), with atomic and molecu- 
lar opacities appropriate to the cooler environments 
of planetary atmospheres and brown dwarfs (Sharp 
& Bm-rows 2007; Burrows & Sharp 1999; Burrows 
et al. 2001, augmented by Burrows et al. 2002 and 
Burrows et al. 2005) . Irradiation in planetary atmo- 
sphere models is incorporated using Kurucz model 
stellar spectra, interpolated to the temperatures 
and surface gravities appropriate to our study (Ku- 
rucz 1979, 1994, 2005). 

In order to examine the effect of metallicity on 

The analysis of Showman et al. 2009 appears to agree 
with this conclusion. 
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planet atmospheres, we compute some models with 
solar abundance of metals, and others with 30 x solar 
abundance. For both metallicities, we derive the 
chemical equilibrium abundances (as functions of 
temperature and density) of several hundred atomic 
and molecular species (Sharp & Burrows 2007). Us- 
ing chemical equilibrium abundances, we calculate 
the monochromatic opacities that are used in the 
radiative transfer calculations. Furthermore, we re- 
compute the equation of state tables for the chemi- 
cal abundances derived for 30 x solar metallicity. 

To treat the redistribution of incident stellar flux 
in a planet's atmosphere, we use the P„ formal- 
ism described in Burrows et al. (2006) and Burrows 
et al. (2008a). In this formalism, in lieu of calcu- 
lating the full three dimensional general circulation, 
the proportion of day-side incident flux that is trans- 
ported to and reradiated from the night-side is pa- 
rameterized as P„, which plausibly ranges between 

0, corresponding to all stellar flux being instanta- 
neously reradiated, and 0.5, corresponding to the 
night side receiving approximately the same amount 
of stellar energy as the dayside (as a result of advec- 
tive heat redistribution). The redistribution takes 
place between 0.003 bars and 0.6 bars. 

For each model planet, we calculate the 

temperature-pressure profile, the planet-star flux 
ratio of emergent infrared radiation, the optical 
albedo, and the temperature and pressure, as func- 
tions of wavelength, of the t\ = 2/3 photosphere. 
It is worthwhile to note that, since photons are not 
stamped as being reflected/ scattered stellar photons 
or emitted photons that arose in a planet's atmo- 
sphere, the albedo is simply the quotient of emer- 
gent flux with incident flux (Burrows et al. 2008b; 
Sudarsky et al. 2003, 2000). Defined this way, it is 
in principle possible for albedos to be greater than 

1. * However, at wavelengths shorter than 1 /zm, the 
albedos we report may be understood as being pre- 
dominantly measures of reflectance and scattering 
properties of the day-side atmospheres. 

3. PLANET MODELS 

The known exoplanets in Neptune's mass range 
are found in a variety of conflgurations. Table 1 
summarizes their diversity. Orbital separations 
range from 0.02 AU to 2.7 AU, and stellar types 
from M through G are represented. The diversity 
of observed planet-star systems motivates our inves- 
tigation of a comparably diverse model set. 

Our goal in this investigation is twofold: («) to 
examine the diagnostics of the known transiting ex- 
trasolar "ice-giants" (GJ436b and HAT-P-llb); and 
{a) to examine more generally the types of radiative 
and thermal structures that might obtain in plausi- 
ble exo-Neptunian atmospheres. To these ends, we 

* In fact, this happens in the mid-infraj:ed, where planets 
are strong thermal emitters. 



compute 29 model atmospheres, listed in Table 2.^'^ 

Six of the models of GJ436b are named 
T35a03Z01P0k0G through T35a03Z01P5k0G (cor- 
responding to a GJ436b-like planet around a GJ436- 
like star, with P„ ranging between 0.0 and 0.5).^ 
Two additional models are T35a03Z01P3k2G (cor- 
responding to GJ436b with P„ = 0.3 and an addi- 
tional optical absorber of opacity k = 0.2 cm^ g~^ 
in the upper atmosphere, where n is the same as Kg 
from Burrows et al. 2008a) and T35a03Z30P0k0G 
(corresponding to GJ436b with P„ = 0.0 and no ex- 
tra optical absorber, and with opacity correspond- 
ing to 30 X solar metallicity). As in previous work, 
K is gray over the range 3 x 10^^-7 x 10'^'* Hz. 
The models of HAT-P-llb are T48a05Z01P0k0H 
through T48a05Z01P5k0H (HAT-P-llb with P„ 
ranging from 0.0 to 0.5) and T48a05Z01P3k2H (ex- 
tra upper atmosphere absorber of opacity k = 
0.2 cm^ g^^)- The additional optical opacity in 
models T35a03Z01P3k2G and T48a05Z01P3k2H is 
conflned to the pressure range from ~5 mbar to the 
top of the atmosphere. 

We also test more generic "ice-giant" models, 
broadly sampling a wide range of possible condi- 
tions. In particular, we test the effect of orbital sep- 
aration, of metallicity, and of stellar type. OGLE- 
05-169L-b, at 2.7 AU, is a bit of an outlier and will 
not be observed in the near future. The remainder 
of the known Neptune-mass exoplanets arc predom- 
inantly at separations a ^ 0.5 AU. This motivates 
our decision to examine planets in this range of sep- 
arations in this initial study. Models testing the 
effects of orbital separation and of metallicity are 
T58a03Z01PlkOG through T58a50Z01PlkOG (solar 
metallicity models from 0.03 AU to 0.50 AU around 
a sun-like star), and T58a03Z30PlkOG through 
T58a50Z30PlkOG (30 x solar metallicity, at the 
same distances from the same star). Models test- 
ing the effect of stellar type are T35a03Z01PlkOG 
through T62a03Z01PlkOG. These models place a 
GJ436b-size planet at 0.03 AU from, respectively, 
GJ436 (T = 3500 K, logioff = 4.80), an MOV star 
(T = 4200 K, logio5 = 4.70), HAT-P-11 (T = 
4780 K, logio5 = 4.60), a KOV star (T = 5200 K, 
logjog = 4.50), a sun-like star (T = 5778 K, 
logiog = 4.44), and an F9V star (T = 6200 K, 
logioff = 4.25). 

^ There are 31 rows in Table 2, of which two, the 26th and 
SOthre duplicates. 

^ These models are publicly available at 
http;//www.astro. princeton.edu/~dsp/cxoneptunes/ and at 
http;/ /www. astro. princeton.edu/~burrows/. 

In this labeling convention, the symbols have the follow- 
ing meanings: "T", for temperature, is followed by 0.01 times 
the star's effective temperature; "a", for semimajor axis, is 
followed by 100 times the orbital semimajor axis in AU; "Z", 
for metallicity, is followed by the ratio with resiject to solar 
abundance of metals; "P" is followed by 10 times P„; "k" is 
followed by 10 times k (in cm^ g^"'); the final letter ("G" or 
"H" ) indicates whether the radius and surface gravity of the 
model planet match those of GJ436b or HAT-P-llb. 
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4. RESULTS 

For each model planet, we present two observ- 
able variables and three that are not observable. 
The observables are the emergent infrared flux, pre- 
sented as a planet-star flux ratio, and the optical 
albedo. The three variables that are unobservable, 
but arc; nonetheless useful probes of our models, 
are the following: the temperature-pressure profile 
(temperature as a function of depth in the atmo- 
sphere); the temperature, as a function of wave- 
length, of the Tx — 2/3 photosphere; and the pres- 
sure, as a function of wavelength, of the same pho- 
tosphere. 

4.1. GJ436b and HAT-P-llb 

The most salient difference between GJ436b 
and HAT-P-llb is that the latter orbits a hotter, 
brighter star. The temperature- pressure profiles of 
models of them, portrayed in Fig. 1, reflect this: 
HAT-P-llb (bottom panel of Fig. 1) is consistently 
hotter at a given pressure. For both profiles, the re- 
distribution of day-side irradiation to the night side 
reduces the dayside's temperature. As P„ increases, 
the thermal profile becomes cooler at a given pres- 
sure in the range where the redistribution takes 
place (0.003 bars to 0.6 bars in these models). This 
effect is more pronounced in GJ436b than in HAT- 
P-llb. 

On the other hand, the effect of an additional 
optical absorber (green curves, k = 0.2 cm^ g~^) is 
more pronounced in HAT-P-llb's profile - precisely 
because its star is hotter and, therefore, brighter in 
the optical. In GJ436b's model, an additional opti- 
cal absorber at high altitude increases upper atmo- 
sphere temperatures by a modest ^60 K. In HAT- 
P-llb's model, however, an extra optical absorber 
increases the temperature at ^^10"^ bars by ~300 K, 
and creates a thermal inversion of about 100 K. 

In the figure for GJ436b, we have included one 

additional model (T35a03Z01P0k0G), with P„ = 
0.0 and 30 x solar opacity. The profile of this model 
looks qualitatively different from those of the other 
models, because it does not quite extend to the 
marginally-stable convection zone. However, the 
Rosseland mean optical depth at the base of this 
model is ~80, and the net flux is constant with alti- 
tude, so the unmodeled portions of the atmosphere 
do not affect the emergent spectrum. 

Let us consider our four wavelength-dependent 
model diagnostics, presented in Fig. 2 (GJ436b) and 
Fig. 3 (HAT-P-llb). We have smoothed the curves 
in each plot to a spectral resolution of 100 for ease 
in viewing. In the flux-ratio plot of Fig. 2 (upper 
left panel), we include the 8-/im measurement that 
Doming et al. (2007a) obtained with the Spitzer In- 
fraRed Array Camera (IRAC). 

The models' emergent infrared fluxes show var- 



ious spectral features. The opacity database con- 
tains significant CH4 cross section at 3.3 /im and 
7.8 /im, H2O and CO absorption at 4.5 /im, and 
H2O at 6 /im and 10 fim (Sharp & Burrows 2007). 
The emergent flux is insensitive to P„ at some 
wavelengths, and shows a moderate-to-strong de- 
pendence on Pn at others (stronger in the models 
of GJ436b than in those of HAT-P-llb). The effect 
of increasing P„ is to reduce the emergent flux in 
"emission" features, speciflcally at ~4 /zm, ~5 /xm, 
~6.5 /xm, and ^10-15 /im. 

The optical albedo (upper right panels of both 
figures) does not exhibit a strong dependence on P„ . 
At wavelengths shorter than 0.5 /<m, the albedo is 
nonnegligible, due to Rayleigh scattering. At longer 
wavelengths, it tends to be ^ 0.05. In both fig- 
ures, the albedos of models with k = 0.0 cm^ g~^ 
start to increase longward of ~0.8 /xm. Albedos of 
HAT-P-llb models are somewhat higher at these 
wavelengths. Note, for both models, the dip in op- 
tical albedo at ~0.95 /xm, which is due to a water 
absorption feature (Sharp & Burrows 2007). 

Although the temperature and pressure of the 
T\ = 2/3 photosphere are not observable (bottom 
left and right panels, respectively), it is instruc- 
tive to examine the effect of increasing P„ on these 
variables. When P„ = 0.0, the temperature varies 
significantly with wavelength, as narrow windows 
of low opacity allow flux from the deeper, warmer 
portions of the atmosphere to escape. As P„ in- 
creases, the location of the photosphere becomes 
more nearly constant with wavelength. The pres- 
sure of the photosphere does not vary much with 
P„, except near 5 /im. As seen in the bottom right 
panel, most of the infrared spectrum originates at 
pressures between 10~^ and 0.1 bars. For these pres- 
sures, as seen in the upper left panel, the day-side 
energy sink for P„ = 0.5 (and to a lesser extent 
for lower P„'s) decreases the T gradient so that for 
P„ = 0.5 there is a quasi-plateau in temperature 
plateau these pressures. 

The effect of an added source of opacity in the 
high atmosphere has striking consequences for both 
the emergent flux and albedo. Although models 
T35a03Z01P0k0G through T35a03Z01P5kOG (k = 
0.0 cm^ g^^) all produce lower planet flux at 
8 /im than Deming et al. (2007a) measured, model 
T35a03Z01P3k2G (k ^ 0.2 cm^ g-i) has an IRAC- 
bandpass-integrated 8-/im flux that is more than 
20% higher than that of T35a03Z01P3kOG (the 
Pn = 0.3 model with no extra optical absorber), and 
is closer to being consistent with the data point's er- 
ror bars.* Furthermore, model T35a03Z30P0k0G, 
with zero redistribution and 30 x solar opacity, has 
an integrated 8-/im flux that is 7.5% higher than 
that of the analogous model with solar opacity. 
Still, its 8-/im flux is 12% lower than bottom of 

* A model with k = 0.5 cm^ (not shown) passes right 
through the center of the error bars. 
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the 1-(J range on the Deming et al. (2007a) mea- 
surement. In both models T35a03Z01P3k2G and 
T48a05Z01P3k2H (GJ436b and HAT-P-llb with 
K = 0.2 cm^ g'^), an extra optical absorber causes 
the upper atmosphere to be warmer and the planet 
flux and the photosphcric temperature to show 
fewer spectra features. Despite the warmer up- 
per atmosphere temperatures and generally warmer 
photospheres, however, the additional optical opac- 
ity causes the the albedo to be much lower in the 
optical ( <i 0.01 throughout much of the optical) be- 
cause Rayleigh scattering is diminished. This is the 
case even though the thermal component of emis- 
sion is greater by dint of the higher temperature.^ 

Consider again the flux-ratio comparisons for 
GJ436b (upper left panels of Fig. 2). All 8 mod- 
els have approximately the same flux at 4.6 /xm, 
but flux varies signiflcantly at 4.0 /xm. The ratio 
of the planet's flux at 4.0 fim to that at 4.6 /im 
varies from -2.4 for T35a03Z01P0k0G to ~1.2 for 
T35a03Z01P4k0G and T35a03Z01P3k2G to ~0.8 
for T35a03Z01P5kOG. Similarly, the ratio of flux 
at ~10 to flux at 4.6 /xm provides a con- 
straint on the redistribution of heat. This ra- 
tio varies from ~7.6 for T35a03Z01P0k0G to ~4.5 
for T35a03Z01P5k0G. These flux ratios, therefore, 
could provide a handle on the amount of redistri- 
bution of heat. In addition, at 7 to 8 /im, extra 
upper atmosphere opacity (T35a03Z01P3k2G) in- 
creases the emergent flux to ~20% above what it 
would be without the extra absorber. We note that 
the difference between the 30 x solar opacity model 
and those with solar opacity appears greater when 
one considers the full spectral shape than when in- 
tegrated over the 8-/um IRAC band. 

The models of HAT-P-llb exhibit a qualitatively 
analogous dependence on P„, although the ratios 
are quantitatively somewhat different because of 
the differences between the two stars and the two 
orbital separations. But, again, T48a05Z01P0k0H 
has significantly greater flux at 4.0 ^m (by a factor 
of -^2) and at 10 yum (by a factor of ~1.4) than 
T48a05Z01P5kOH. Furthermore, because HAT-P- 
llb is a hotter star with correspondingly more of 
its flux in the optical, an extra optical absorber in 
the upper atmosphere has a more significant effect 
on this planet than on GJ436b. Between 7-8 /xm, 
T48a05Z01P3k2H is approximately twice as bright 
as the models without an extra absorber, and is also 
much brighter (by ~40%) at wavelengths greater 
than 15 /xm. 



4.2. Distance 

The temperature of an irradiated stellar com- 
panion's photosphere drops off rapidly with increas- 

^ At sufficiently liigli values of K and with sufficient irra- 
diation, this trend would reverse, and higher k would lead to 
a larger apparent albedo. 



ing orbital separation, and its emergent fiux de- 
creases correspondingly. For separations greater 
than ~0.5 AU, both the reflected/scattered optical 
light and the infrared emission become too dim to be 
detected without spatially separating (and masking) 
the light of the star. Although direct detection of 
wide-separation Neptune-mass objects will be feasi- 
ble in the future (Burrows 2005), in this paper we 
focus on objects within 0.5 AU of their primaries. 

The top panel of Fig. 4 presents temperature- 
pressure profiles of models at 0.03, 0.05, 0.10, 
0.30, and 0.50 AU from a sun-like star (G2V). 
In these models (T58a03Z01P0k0G through 
T58a50Z30P0k0G in Table 2), k = 0.0 cm^ g-i 
and Pn = 0.0. Solar opacity models are represented 
by solid curves, and 30 x solar opacity models with 
dashed curves; the differences between these are 
discussed in §4.3 below. 

Figure 5 contains four panels analogous to the 
four of Fig. 3. The predominant trends with dis- 
tance in the flux-ratio plot (upper left panel) and in 
the photosphcric temperature plot (lower left panel) 
are, respectively, that both the planet fiux and 
planet temperature decrease as the model planet is 
moved farther from the star. In addition, there is 
some change in the spectral shape - cooler, more 
distant planets have redder spectra. The albedos of 
these models (upper right panel of Fig. 5) change 
from low at short wavelengths and high (—0.7) at 
long wavelengths for the 0.03 AU models to mod- 
erate (~0.3) at short wavelengths and low at long 
wavelengths for the 0.50 AU models. The close- 
in models have relatively high apparent albedo at 
wavelengths approaching 1 /xm because of the con- 
tribution from the Wien tails of their thermal emis- 
sion. 



4.3. Metallicity 

It was initially surprising to us that models 
with higher metallicity around solar-type stars, with 
their concomitant changes in both opacity and equa- 
tion of state, nevertheless have extremely similar 
spectral profiles to their lower metallicity cousins 
(solid and dashed lines in top left panel of Fig. 5). 
The temperature-pressure profiles (top panel of 
Fig. 4) of 30 X solar opacity models are quite dif- 
ferent from those of solar opacity models; increas- 
ing metallicity appears to translate the T P profiles 
up to higher altitudes (lower pressures). But the 
photosphcric temperature (bottom panel of Fig. 5) 
remains essentially unchanged with higher opacity. 

How is this possible? The photosphere pressure 
plot (lower right panel of Fig. 5) offers the explana- 
tion. Because, viewed from the outside, the optical 
depth increases so much more quickly with depth 
into atmosphere for the high metallicity models, 
their tx = 2/3 surfaces are at significantly lower 
pressures (by roughly an order of magnitude), cor- 
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responding to much higher altitudes. But, as Fig. 4 
shows, higher in the atmosphere the temperature 
is higher in the high metallicity models, and these 
two effects conspire to leave the temperature of the 
photospheres of the high metallicity models virtu- 
ally unchanged (bottom left panel of Fig. 5). 

The changes to the albedos (top right panels of 
Fig. 5) are not dramatic, although at the short wave- 
length end the albedos of the 30 x solar metallicity 
models are about half that of the solar metallic- 
ity models. Furthermore, at the long wavelength 
end, a deep absorption feature due to water arises 
at ~0.95 /Ltm in the high metallicity models that is 
not present in the low metallicity ones. Since the 
two characteristics that are significantly different 
at high metallicity (the T-P profile and the photo- 
spheric pressure) are not observable, these compara- 
tively modest changes to albedo might represent our 
best chance to constrain the atmospheric metallic- 
ity through whole-disk observations (though transit 
spectrum observations could provide a complemen- 
tary constraint). Note, however, that this water fea- 
ture is the same one seen in solar metallicity models 
of GJ436b and HAT-P-llb in Figs. 2 and 3, and, 
therefore, might be a signature of high metallicity 
only in the context of a modeling effort that includes 
the appropriate stellar irradiance spectrum. 

For a solar-type irradiation spectrum, the opac- 
ity (solar or 30 x solar) of model planet's atmo- 
sphere has a minor influence on the emergent spec- 
trum. However, the opacity has a stronger influ- 
ence when the irradiation spectrum is from a cool 
M-dwarf (as shown in Fig. 2, where the irradiation 
is from GJ436). 

4.4. Stellar Type 

How do a Neptune- mass planet's structure and 
appearance depend on its host star? Models 
T35a03Z01P0k0G through T62a03Z01P0k0G ad- 
dress this question. These models are of planets the 
size of GJ436b. separated 0.029 AU from stars rang- 
ing from cool (3500 K) to moderately hot (6200 K). 

The most prominent differences between models 
around cooler and dimmer stars and those around 
hotter stars is that, unsurprisingly, planets around 
hotter stars are hotter and bluer. The bottom panel 
of Fig. 4 and the bottom left panel of Fig. 6 indi- 
cate how more massive, brighter stars lead to hotter 
planets at a given distance. 

For a planetary system that is of order a Gyr 
or older, a Neptune-mass planet has lost almost 
all of its heat of formation (Burrows et al. 2000, 
2003; Baraffe et al. 2005). As a result, for planets 
at separations <0.5 AU, contributions to the emer- 
gent planetary flux (reflected/scattered and emit- 
ted) from self-luminance are miniscule. Therefore, 
the integrated planet flux divided by the integrated 
stellar flux is, to a high degree of approximation. 



simply 1 / (47r) times the solid angle subtended by 
the planet from the vantage of the star (np/47r = 
wRp^ /Atto^). This leads to the somewhat counter- 
intuitive flux-ratios plot in the upper left corner of 
Fig. 6: in the mid-infrared (at wavelengths greater 
than ^8 /im), a planet around a cool star is signifi- 
cantly brighter relative to its star than one around 
a hot star. Of course, this is because the same total 
relative flux is distributed predominantly at longer 
wavelengths for a cooler planet and at shorter wave- 
lengths for a hotter planet. 

Interestingly, at wavelengths less than ^-^8 /im, 
the photospheric pressures are not dramatically dif- 
ferent for the six models, dcispitc very different T- 
P profiles. Finally, similar to the trend seen in 
albedo versus distance seen in §4.2, cooler models 
have higher albedos at short wavelengths, while hot- 
ter models have higher apparent albedos at longer 
wavelengths as the Wien tail enters into the optical 
and near infrared. 

5. SUMMARY AND CONCLUSIONS 

We have presented a series of theoretical 1- 
D models of the radiative and thermal structures 
of Neptune-mass, "ice-giant" exoplanets, including 
models at opacities corresponding both to solar 
metallicity and to 30 x solar. To produce the latter, 
we computed equilibrium chemical abundances of 
hundreds of molecular and atomic species, and we 
recomputed the equation of state with the higher 
mean molecular weight. We produced models of 
the transiting planets GJ436b and HAT-P-llb, and 
models of generic planets in a range of plausible 
conditions. 

In our investigation of GJ436b and HAT-P-llb, 
we find the following: 

• The ratio of planet flux at 4.0 /im or at 10 /im 
to that at 4.6 /im might place a constraint 
on the amount of redistribution of heat from 
the day to the night side: more redistribu- 
tion leads to a smaller flux ratio. The ratio of 
planet flux at 7 /im to that at 4.6 /im might 
place a constraint on the possible presence of 
an extra optical absorber in the high atmo- 
sphere. 

• An extra optical absorber, if present in either 
planet, would have a more dramatic effect on 
both the T-P profile and the emergent spec- 
trum of HAT-P-llb than of GJ436b, because 
HAT-P-11 is brighter in the optical. 

• Models of GJ436b without an extra ab- 
sorber are significantly dimmer than the Dem- 
ing et al. (2007a) measurement, while the 
model with an extra absorber of opacity 
K = 0.2 cm^ g~^ is closer to being consis- 
tent with the measurement, though still shy 
of it. 
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• Models of these two planets have moderate 
optical albedos (~0.3) at the short wavelength 
end of the spectrum, and low albedos at longer 
wavelengths. 

Furthermore, we find the following broad trends 
among generic models of planets in this mass-range: 

• Higher metallicity implies higher opacity for a 
given column mass of atmosphere. As a result, 
photospheres of higher metallicity models are 
at lower pressures. 

• Nonetheless, when the irradiation spectrum is 
from a solar-type star, higher metallicity mod- 
els appear to differ little from solar metallic- 
ity models in observable variables. The most 
noticeable difference, among close-in planets 
around solar-type stars, is the water absorp- 
tion feature at ~0.95 /xm that shows up in 
higher metallicity models. When the irradia- 
tion spectrum is from a cool M-dwarf (GJ436), 
the metallicity of the planet's atmosphere has 
a larger influence on the emergent spectrum. 

• Neptune-mass planets around cooler stars 
have larger flux ratios in the mid-infrared than 
ones at the same distance around hotter stars. 



(This conclusion is germane to Jupiter-mass 
planets, as well.) 

Close-in Neptune-mass planets present an ob- 
servational challenge because, with radii a factor 
of ~ 3 smaller than Jupiter-mass planets, they are 
roughly an order of magnitude dimmer. However, 
observations of GJ436b have already been made 
with Spitzer, and additional ones are planned. For 
GJ436b, 3.6 and 4.5 /xm IRAC observations will 
help to quantify both the amount of redistribution 
of heat from the day to the night sides and the pos- 
sible presence of an extra absorber in the high atmo- 
sphere. Future observations with new facilities will 
provide constraints on models of other planets in 
this mass range, and will sharpen our understand- 
ing of this class of objects. 
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TABLE 1 

Characteristics of -^Neptune-Mass Exoplanets 



Planet 


Mass 


Radius 


logio 9 


Semi-major Axis 


Stellar Type 


References 




(J^Nep'') 


(cgs) 


(AU) 


CoRoT-7b 


0.65 


0.43 


3.59 


0.02 


KO V 


Pressin et al. 2009 


GJ436b 


1.33 


1.04 


3.15 


0.03 


M2.5 


Torres 2007 


HAT-P-llb 


1.50 


3.05 


1.22 


0.05 


K4 


Bakos ct al. 2009 
Dittman ct al. 2009 


HD 181433b 


0.44 


- 


- 


0.08 


K3 IV 


Bouchy et al. 2009 


HD 285968b 


0.49 


- 


- 


0.07 


M2.5 V 


Forveille et al. 2009 


HD 40307d 


0.53 


- 


- 


0.13 


K2.5 V 


Mayor et al. 2009 


HD 7924b 


0.54 


- 


- 


0.06 


KO V 


Howard et al. 2009 


HD 69830b 


0.61 


- 


- 


0.08 


KO V 


Lovis et al. 2006 


HD 160691c 


0.62 


- 


- 


0.09 


G3 IV-V 


Pcpc ct al. 2007 


55 Cnc e 


0.63 


- 


- 


0.04 


G8 V 


Povcda & Lara 2008 


GJ 674b 


0.69 


- 


- 


0.04 


M2.5 


Bonfils et al. 2007 


HD 69830c 


0.70 


- 


- 


0.19 


KO V 


Lovis et al. 2006 


OGLE-05-169L b 


0.74 


- 


- 


~2.7 


- 


Gould et al. 2006 


HD 4308b 


0.87 


- 


- 


0.11 


G5 V 


O'Toole et al. 2009b 


GJ 581b 


0.91 


- 


- 


0.04 


M3 


Mayor et al. 2009 


HD 190360c 


1.06 


- 


- 


0.13 


G6 IV 


Vogt et al. 2005 


HD DaooUd 


1 no 
i.Uo 






U.DO 


-K.U V 


Lovis ct al. 2UUD 


HD 219828b 


1.22 






0.05 


GO IV 


Melo et al. 2007 


HD 16417b 


1.28 






0.14 


Gl V 


O'Toole et al. 2009a 


HD 47186b 


1.33 






0.05 


G5 V 


Bouchy et al. 2009 


HD 99492b 


2.02 






0.12 


K2 V 


Marcy et al. 2005 


HD 11964b 


2.04 






0.23 


G5 


Wright ct al. 2009 


HD 49674b 


2.13 






0.06 


G5 V 


Wright ct al. 2007 


55 Cnc f 


2.67 






0.78 


G8 V 


Fischer et al. 2008 



s-The mass of Neptune is MNep = 1.02 x lO^^ g = 0.054 Mjup, where Mjup = 1.90 x 10^3 g. 
''The radius of Neptune is ilNep = 2.48 x 10^ cm = 0.346 i?jup, where Rj^p = 7.15 x 10^ cm. 
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TABLE 2 
ExoPLANET Models* 



Model Name 


Stellar Model 








a 

(AU) 


Mctallicity 
(Zq) 


Pn 


Rp 

(-Rncp) 


logio 9 
(cgs) 


K 

(cm^ g-i) 


T35a03Z01P0k0G 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T35a03Z01PlkOG 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.1 


1.04 


3.145 


0.0 


T35a03Z01P2kOG 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.2 


1.04 


3.145 


0.0 


T35a03Z01P3kOG 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.3 


1.04 


3.145 


0.0 


T35a03Z01P4kOG 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.4 


1.04 


3.145 


0.0 


T35a03Z01P5kOG 


GJ436 (T = 3500 K 


log 


10 9 


= 4.80) 


0.03 


1 


0.5 


1.04 


3.145 


0.0 


T48a05Z01P0k0H 


HAT-P-11 (T = 4780 K; log 


10(7 = 4.60) 


0.05 


1 


0.0 


1.22 


3.052 


0.0 


T48a05Z01PlkOH 


HAT-P-11 (T = 4780 K 


log 


10 5 = 4.60) 


0.05 


1 


0.1 


1.22 


3.052 


0.0 


T48a05Z01P2kOH 


HAT-P-11 (T = 4780 K 


logio 9 = 4.60) 


0.05 


1 


0.2 


1.22 


3.052 


0.0 


T48a05Z01P3kOH 


HAT-P-11 (T = 4780 K 


log 


10 5 = 4.60) 


0.05 


1 


0.3 


1.22 


3.052 


0.0 


T48a05Z01P4kOH 


HAT-P-11 (T = 4780 K 


log 


10 5 = 4.60) 


0.05 


1 


0.4 


1.22 


3.052 


0.0 


T48a05Z01P5kOH 


HAT-P-11 (T = 4780 K 


log 


10 5 = 4.60) 


0.05 


1 


0.5 


1.22 


3.052 


0.0 


T35a03Z01P3k2G 


GJ436 (T = 3500 K; log 


10 S 


= 4.80) 


0.03 


1 


0.3 


1.04 


3.145 


0.2 


TlSaOr)Z01P3k2II 


HAT-P-11 (T = 1780 K 


log 


i„ '/ = 1-00) 


0.05 


1 


0..3 


1.22 


zm2 


0.2 


T35a03Z30P0k0G 


GJ436 (T = 3500 K; logio 9 


= 4.80) 


0.03 


30 


0.0 


1.04 


3.145 


0.0 


T58a03Z01P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T58a05Z01P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.05 


1 


0.0 


1.04 


3.145 


0.0 


T58al0Z01P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.10 


1 


0.0 


1.04 


3.145 


0.0 


T58a25Z01P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.25 


1 


0.0 


1.04 


3.145 


0.0 


T58a50Z01P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.50 


1 


0.0 


1.04 


3.145 


0.0 


T58a03Z30P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.03 


30 


0.0 


1.04 


3.145 


0.0 


T58a05Z30P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.05 


30 


0.0 


1.04 


3.145 


0.0 


T58al0Z30P0k0G 


G2 V (T = 5778 K; 


logic 9 = 


= 4.44) 


0.10 


30 


0.0 


1.04 


3.145 


0.0 


T58a25Z30P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.25 


30 


0.0 


1.04 


3.145 


0.0 


T58a50Z30P0k0G 


G2 V (T = 5778 K; 


logio 9 = 


= 4.44) 


0.50 


30 


0.0 


1.04 


3.145 


0.0 


T35a03Z01P0k0G'' 


GJ436 (T = 3500 K; log 


10 


= 4.80) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T42a03Z01P0k0G 


MO V (T = 4200 K; logio 9 = 


= 4.70) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T48a03Z01P0k0G 


HAT-P-11 (T = 4780 K; log 


10 5 = 4.60) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T52a03Z01P0k0G 


KO V (T = 5200 K; logy, g = 


= 4.50) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T58a03Z01P0k0G'= 


G2 V (T = 5778 K; logio 9 = 


= 4.44) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 


T62a03Z01P0k0G 


F9 V (T = 6200 K; logio 9 = 


-- 4.25) 


0.03 


1 


0.0 


1.04 


3.145 


0.0 



This table contains 31 rows describing the 29 models we computed. Two rows are repeated. 

*The model name is a shorthand for the model characteristics, as follows: "T" is for temperature, followed by the temperature 
in Kelvin divided by 100; "a" is for semimajor axis, followed by 10 times the orbital separation in AU; "Z" is for metallicity, 
followed by the multiple of solar abundance; "P" is for P„, followed by 10 x P„; "k" is for k, followed by 10 x k (in cm^ g~^). 
The final letter is either "G" or "H" : if "G" , the planet's radius and surface gravity are those of GJ436b; if "H" , they are those 
of HAT-P-llb. 

''This is a repeat of the first row. 

''This is a repeat of the 15th row. 
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600 800 1000 1200 1400 1600 1800 2000 2200 

T(K) 



Fig. 1. — Model temperature-pressure profiles of GJ436b and HAT-P-11 b. Top: Models of GJ436b, with the redistribution 
parameter P„ varying from (no redistribution) to 0.5 (complete redistribution of heat to the night side). A model is shown 
with added upper atmosphere opacity k = 0.2 cm^ g~^ (here, k is the same as the Ke from Burrows et al. 2008a). An additional 
model shows the influence of SOxsolar atmospheric opacity. Bottom: Same, for HAT-P-11 b, except without the SOxsolar 
model. 
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Planet-Star Flux Ratios 
GJ436b Models 



P = 

n 

P = 

n 

P = 

n 

P = 

n 

P = 

n 

P = 

n 

P = 



0.0; 30x solar 
0.5 Deming et al. (2007) 
0.4 
0.3 
0.2 
0.1 
0.0 

0.3; K = 0.2 




Temperature of x =2/3 Photosphere 



GJ436b Models 




0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



P = 0.0; 30x solar P =0.3;k = 0.2 

n ,1 



10"' 



10 



I (nm) 



10' 




Pressure of x^=2/3 Photosphere 
GJ436b Models 



10 



X (|jm) 



Fig. 2. — Additional characteristics of models of GJ436b depicted in the top panel of Fig. 1. Top left: Planet-star flux 
ratios. Varying P„ has little influence on the planet's flux at 8 fim, but an additional upper atmosphere absorber results in a 
significantly better match of the model to the data. 30 X solar atmospheric opacity also helps bridge the gap between the solar- 
opacity models and the observed data. The black square is the Spitzer IRAC-4 measurement of Deming et al. (2007a); vertical 
bars indicate l-cr uncertainty; the horizontal bar indicates the full width at 10% maximum of the IRAC-4 band's transmission 
function. Large X's indicate integrated planet flux divided by integrated stellar flux over this band. Top right: Albedos. The 
albedo is lower for lower P„, and much lower for nonzero k. Bottom left: Temperature of the tx = 2/3 surface, as a function 
of wavelength. Bottom right: Pressure of the tx = 2/3 surface, as a function of wavelength. 
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Planet-Star Flux Ratios 
HAT-P-1 1b Models 




Temperature ofT^=2/3 Photosphere 
HAT-P-1 1b Models 

P : 
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^ = 0.1 
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P 


' = 0.0 


Ah I V V 1 /iir 
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Pressure of x^=2/3 Photosphere 






HAT-P-1 lb Models 



10 



I (nm) 
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Fig. 3. — Additional characteristics of models of HAT-P-llb depicted in the bottom panel of Fig. 1. Completely analogous 
to Fig. 2. 
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T-P Profiles, "Neptune" Models 

Solar opacity: solid 

solar opacity: dashed 
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Fig. 4. — Model temperature-pressure profiles of GJ436b-like planet. In all models, the redistribution parameter P„ = 0. 
Top: Models are computed at a range of distances from a Sun-like star, at solar and at 30-times solar metallicity. More distant 
models are cooler, and more metal-rich models are warmer at a given pressure. Bottom: Models are computed for a range of 
host stars, from GJ436 (a cool M-dwarf) through an F9 star (6200 K). 
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Fig. 5. — Additional ciiaracteristics of models for different orbital distances and metallicities depicted in the top panel of 
Fig. 4. Analogous to Fig. 2. Top left: Planet-star flux ratios for a GJ436b-like model planet at a variety of distances from 
a G2V star, at solar and at 30-times solar metallicity. Changes in metallicity have only a very modest influence on model 
spectra. Top right: Optical albedos for these same models. At high metallicity, generally the cooler, more distant models 
(dashed curves) have larger optical albedo. At solar metallicity, the albedo varies less smoothly with semimajor axis. Bottom 
left: Temperature of the t\ = 2/3 surface, as a function of wavelength. Bottom right: Pressure of the = 2/3 surface, as a 
function of wavelength. 
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Fig. 6. — Additional characteristics of models for different host stars 
Fig. 2. 
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depicted in the bottom panel of Fig. 4. Analogous to 



